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Abstract 

Several distinct colour flow topologies are possible in multiparton config- 
urations. A method is proposed to find the correct topology, based on a 
minimization of the total transverse momentum of produced particles. This 
method is studied for three-jet Z° — > qqg and four-jet W + W~ — ► qiq^zQ-A 
events. It is shown how the basic picture is smeared, especially by parton- 
shower activity. The method therefore may not be sufficient on its own, but 
could still be a useful complement to others, and e.g. help provide some 
handle on colour rearrangement effects. 
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When high-energy processes produce multiparton states, it is generally believed that 
the confinement property of QCD leads to the formation of colour flux tubes or vortex lines 
spanned between the partons. These tubes/vortices are here called strings, in anticipation 
of our use of hadronization based on the string model A quark or antiquark is attached 
to one end of a string. A gluon is attached to two string pieces, one for its colour and one 
for its anticolour index, and thus corresponds to a kink on the string. 

The simplest kind of events, e + e~ — ► "y*/Z° — > qq, gives a single string stretched be- 
tween the q and the q. The direction of the colour flow can, in principle, be distinguished 
by flavour correlations ||, but will not be studied here. In next order, qqg events corre- 
spond to a string stretched from the q via the g to the q. The colour topology is unique, 
but experimentally it is not normally known which of the three jets is the gluon one, so 
this gives a threefold experimental ambiguity. 

From four partons onwards, true ambiguities of the topology exist, even when the 
identity of the partons is known. In qqgg events, the string can be drawn from the q to 
either of the two gluons, on to the other gluon and then to the q. This gives two possible 
topologies. A third topology, not expected to leading order in N c , is when one string runs 
directly between the q and q and another string in a closed loop between the two gluons. 
There is an experimental ambiguity, in picking the two quarks among the four partons, 
which gives a further factor of six, i.e. a total of eighteen possible topologies (reduced to 
fifteen if single and double strings are not distinguished). 

Another four-jet final state is obtained in the process e + e~ — ► W + W~ q{q 2 ( lzQ.^ 
Since the W's are colour singlets, in principle each of qiq 2 and q^q^ form a separate colour 
singlet string. However, by soft gluon exchange or some other mechanism, alternatively 
qiq A and q^q 2 may form two singlets. Since it would not normally be known which of the 
four jets are quarks and which antiquarks, there is a total of three experimental pairings 
of four jets, where the third corresponds to the unphysical flavour sets q\qz and q 2 q 4 . 

Among the topologies above, only the three-jet qqg events have been studied in detail. 
It has been shown that the string approach here correctly predicts the topology of particle 
flow, with a dip in the angular region between the q and q jets ||, |J. This comes about 
because the qg and qg string pieces produce (soft) particles in the respective angular 
ranges, while there is no string directly between the q and q. The same effect is also 
obtained as a consequence of colour coherence in perturbative soft-gluon emission M 
normally these two approaches give the same qualitative picture. Although by now LEP 1 
has produced large samples of four-jet events, the energy flow between these jets have not 
been studied in detail, presumably because of the large number of possible topologies. 

LEP 2 will provide four-jets from W pairs, and here the issue of colour topology may 
become of great importance |6|, |7[]. If, by colour rearrangement, an original q{q 2 plus q 3 q 4 
colour singlet configuration is turned into a q{q A plus qzq 2 one, particle production will be 
somewhat different. Methods to determine the W mass from LEP 2 four-jet events will 
then give different results. There is more than one model of the colour rearrangement pro- 
cess; therefore the uncertainty on the W mass could be as high as 100 MeV M, i.e. larger 
than the expected statistical error of the order of 40 MeV. The final number will then 
be dominated by the mixed hadronic-leptonic channel W + W~ — ► tviqq^ which has about 
the same statistics. The hadronic events could be recuperated if colour rearrangement 
effects could be diagnosed from the data itself. Additionally, the ability to distinguish 
between various reconnection scenarios could provide information on the nature of the 
QCD vacuum, and therefore be of fundamental interest. Unfortunately, realistic recon- 
nection scenarios give only very minute effects on the data (remember that the effect on 
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the W mass is at most of the order of one per mille), and so attempts to find useful signals 
have been rather unsuccessful 0. There is only one claim for having found a signal |J 
where, in a few events, a central rapidity gap separates the particle production from two 
low-mass colour singlet reconnected systems. It has turned out, however, that neglected 
angular correlations in the W pair decays tend to reduce this signal to the border of 



observability [[U], |TT[ . Therefore one would like to devise alternative methods to diagnose 
the appearance of colour rearrangement. 

We now want to propose and study a method to select the correct string topology. The 
starting point is the observation that, were it not for various smearing effects, hadrons 
would be perfectly lined up with the string pieces spanned between the partons. In mo- 
mentum space the hadrons would appear along hyperbolae with the respective endpoint 
parton directions as asymptotes. In a frame where a string piece has no transverse mo- 
tion, i.e. where the two endpoint partons are moving apart back-to-back, the hadrons 
produced by this piece would have vanishing transverse momentum. A hadron could then 
successively be boosted to the rest frame of any parton pair until the pair is found for 
which the particle p± is vanishing. When smearing is introduced, this p± will no longer be 
vanishing, but it still would be reasonable that a "best bet" is to assign each hadron to the 
string piece with respect to which it has smallest p±. The most likely string configuration 
would be the one where the sum of all hadron p±'s is minimal. While the fluctuations for 
each single hadron is too large for usefulness, it could be hoped that the net effect of all 
the hadrons would be to single out the correct topology. 

To be more precise, here is the scheme proposed, to be carried out for each event: 

1. Use some jet clustering algorithm to identify the directions of the number of jets 
that should be used for the current application. 

2. Enumerate the possible colour flow topologies allowed in the process. 

3. Form one J2p± measure for each topology, where the sum runs over all final-state 
particles in the event. For each particle, its p± is defined as the minimum of the 
p±'s obtained by boosting the particle to the rest frame of each of the string pieces 
making up the current topology. 

4. Identify the correct topology as the one with smallest Y^P±- 

Several effects could smear the simple picture and lead to incorrect conclusions. The 
main ones are: 

• Errors in the reconstruction of jet directions. 

• Additional p± caused by perturbative QCD branchings, predominantly gluon emis- 
sion. 

• The p± generated by the fragmentation process. 

• Secondary decays of unstable hadrons. 

A convenient test bed for the relative importance of these effects is provided by qqg 
events. In symmetric three-jet events at Z° energies, without any parton-shower activity, 
the gluon is correctly identified in 87% of the cases. This should be compared with the 
33% expected from a random picking among the three alternatives. When complete Z° 
events events are generated (with Pythia 5.7 and Jetset 7.4 Hl2"|), three jets are found 
(in 10-20% of the events) and the method above is used, the success rate is around 
55-60%. The "correct" answer is here found by tracing the original quark and antiquark 
through the shower history and associating them with the jets they are closest to in angle. 
The conclusion of this kind of studies is that the method does work, though maybe not 
as well as one might have hoped, and that the main cause of errors is the perturbative 
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gluon emission. 

The studies can be extended to four-jet qqgg events, although the large number of 
colour flow topologies here makes the method rather inefficient for selecting the correct 
topology. A more modest objective is to establish differences between string and inde- 
pendent fragmentation ]I3| models. In the latter approach, particle production is aligned 
along the direction connecting a parton with the origin in the cm. frame of the event, 
again with some smearing effects. This approach is already disfavoured by the three-jet 
studies but is a convenient reference. We have not completed a full study, but can il- 
lustrate with results for a simple four-jet cross topology with q and q back-to-back. Fig. |l]a 
gives the difference in J^Pi. between the worst topology (where none of the string pieces 
run the way assumed) and the right one. Results are for charged particles, and the J2p± 
has been normalized to the number iV of charged particles per event. The independent 
fragmentation distribution is symmetric around the origin, as it should, while the string 
approach shows a clear offset. Without knowledge of the correct answer, one is reduced 
to studying a measure such as the difference between the maximal and the minimal P± 
among all the twelve possible topologies, Fig. ||b. By construction this is a positive num- 
ber, also for independent fragmentation, but an additional shift is visible for the string 
approach. Identification of one or both quarks, e.g. by b quark tagging or energy ordering 
(quarks are likely to have higher energy than gluons), would cut down on the number of 
topologies to be compared and therefore enhance the signal. Some experimental studies 
along these lines could therefore be interesting. 

We now turn to the main application of this letter, namely four-jet W + W~ events. 
Since the two W decay vertices are less than or of the order of 0.1 fm apart, while typical 
hadronic distances are of the order of 1 fm, the two W decays occur almost on top of 
each other. QCD interconnection effects could appear at all stages of the process, namely 
in the original perturbative parton cascades, in the subsequent soft hadronization stage, 
and in the final hadronic state. It can be shown that perturbative effects are strongly 
suppressed |7j, but no similar arguments hold for the other two. Bose-Einstein effects in 



the hadronic state could be the largest individual source of W mass uncertainty [14], but 



it is the least well studied. The presence of such Bose-Einstein effects presumably could 
be established from the data itself, while reconnection in the hadronization stage is less 
easy to diagnose. We will study whether the YjPi. measure offers any help here. 
The reconnection models used as references in this work are: 

1. Reconnection after the perturbative shower stage but before the hadronization, with 
reconnection occuring at the 'origin' of the showering systems. This 'intermediate' 
model is the simplest of the more realistic ones. 

2. Reconnection when strings overlap based on cylindrical geometry. A 'bag model' 
based on a type I superconductor analogy. The reconnection probability is propor- 
tional to the overlap integral between the field strengths, with each field having a 
Gaussian fall-off in the transverse direction, the radius being about 0.5 fm. The 
model contains a free strength parameter that can be modified to give any recon- 
nection probability. 

3. Reconnection when strings cross. In this model the strings mimic the behaviour of 
the vortex lines in a type II superconductor, where all topological information is 
given by a one-dimensional region at the core of the string. 

4. Reconnections occur in such a way that the string 'length' is minimized. As a 
measure of this length the so-called A-measure is used, which essentially represents 
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the rapidity range for particle production counted along the string. This can also 

be seen as a measure of the potential energy of the string. 
Models 1 through 3 is described in and the last in ||. Further models have been 
proposed [ITS . 



For the study, events should have a clear four-jet structure. To achieve this we demand 
that each jet must have some minimum energy and that the angle between any jet pair 
must not be too small J7|. When applied to the expected statistics of LEP 2, the number 
of events left after the cuts will be about 2500 per experiment; therefore statistics will be 
a problem when different models are compared with each other. 

Three different algorithms are used to identify which jet pairs belong together: 

1. The qi^QsQ^ configuration before parton showers can be matched one-to-one with 
the reconstructed jets after hadronization. This is done by minimizing the products 
of the four (jet+g) invariant masses. The original quark information is not available 
in an experimental situation, so this measure can only be used as a theory reference. 

2. The invariant mass of jet pairs from the same W should be close to the known W- 
mass of about 80 GeV. Among the three possible jet pairings, therefore the one is 
selected which has minimal \rriij — 80| + |m&j — 80|, where i,j, k, I are the four jets. We 
have picked this method rather than a few similar ones since it has (marginally) the 
best correlation with the reference method above. This method is mainly probing 
the electroweak aspect of W pair production, namely the W mass spectrum, while 
it should be less dependent on the QCD stages of showering and hadronization. 

3. The J2p± method introduced in this letter provides an alternative measure, that 
rather should be sensitive to the QCD stages and less so to the electroweak one. 
Without colour reconnection it should (hopefully) agree with the previous two, while 
it could give interesting differences if reconnection occurs. 

The agreement between these three methods is shown in Table [1], with and without re- 
connection, the former for model 1. As should be expected, algorithm 2 comes close to the 
"correct" answer of number 1, and is not significantly affected by colour rearrangement. 
The J2p± method, algorithm 3, shows the expected dependence on colour rearrangement, 
with a smaller success rate when colour rearrangement is included. Note, however, that 
the success rate does not drop below the naive 33% number, indicating that the J2p± 
method is also picking up other aspects of events, such as the jet topology. The results in 
the table are for a 170 GeV energy, but we do not expect a significant energy dependence. 

Both methods 2 and 3 can be applied to data, so therefore the correlation is an 
observable. The rate of reconnection could be extracted, by interpolation between the 
two extremes of no and complete colour rearrangement. Statistically it should be feasible 
to establish a signal for reconnections, if they occur at a rate above the 10-20% level. 
The systematic errors on the correlation method may be large, however, especially for 
the model-dependent change when reconnection is included. It is therefore important to 
study whether a differential distribution would better highlight qualitative differences. 

Algorithm 2 can be used to identify the best hypothesis for which jets should be 
paired to form the two W's, and also the worst hypothesis, where \rriij — 80| + |rajy — 80 1 
is maximal. The J2p± can be calculated for each of these two extremes. When the 
strings are reconnected, the first sum should increase and the second one decrease relative 
to the no-reconnection case. The signal is therefore enhanced by making use of the 
difference, A = (^P-iJworst — (J2p±)bcst, which should decrease in case of reconnection. 
The subtraction furthermore has the advantage of removing some spurious fluctuations 
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from the comparison. The main example is high-momentum particles, where the assumed 
string hyperbolae attach well with the four jet directions and therefore all three string 
hypotheses give the same contribution to J2P±- (Our studies show that particles with 
momenta above 3 GeV add little to the discrimination between the string hypotheses.) 

The A measure is plotted for models 1-4 in Fig. |^. Note that the results for models 
1 and 4 correspond to 100% reconnection, while the reconnection rate in models 2 and 3 
is about 30%. (The reconnection fraction could be varied in all models, but the effects of 
reconnection are not linear in this fraction for models 2 and 3, so a reasonable value is 
preferred here.) All reconnection models show the expected shift towards smaller A values, 
and the magnitude of the shift is comparable once differences in assumed reconnection 
fractions are removed. Remaining differences imply that one cannot model-independently 
extract a reconnection rate from the data. The signal for reconnection may be enhanced, 
compared with the results of Table [I], by cuts on A, e.g. by only considering the fraction 
of events with A < 0. This is at the price of a reduced statistics, however, so the balance is 
not so clear. It may be better to use measures that gauge the full shape of the curve, given 
that the physics of the no-reconnection scenario is presumed well-known (by extrapolation 
from the Z° results). 

Prospects look promising to diagnose colour rearrangement along these lines but, as 
before, the combination of low effects and small statistics could give marginal results. 
Furthermore, hopes should not be raised too high that this would immediately imply a 
scheme to correct a W mass measurements for the reconnection effects: of the models 
above, number 3 shifts the W mass downwards while the others shift it upwards ]7|, |j. 
and yet they all shift the A distribution in the same direction. Clearly the study of 
reconnection effects ultimately must be based on a host of different measures, the J2p±. 
one and others. 

It may be of some interest to understand why effects are not larger. Several simulations 



have been performed with various simplified toy models to study this issue [lq|. It turns 



out that there are two main mechanisms that smear distributions and make them less 
easily distinguished. One is parton showers, just as for the Z° — > qqg process studied 
above. The other is the geometry of the process, namely that the helicity structure of 
the W + W~ — > ?ig 2 ?3?4 process is such that q\ and q 3 tend to go in the same general 
direction, as do q 2 an d Qa [03- The overall change of string topology by reconnection 
therefore is not as drastic as if the q± had tended to be close to g 4 and vice versa. Had it 
been possible to remove these effects, i.e. study events without shower activity and with q\ 
and q 4 reasonably close in angle, the original and the colour-reconnected A distributions 
would almost completely separate. In practice, only a modest reduction of shower activity 
could be obtained by requiring that all four jets be reasonably narrow, while tagging of 
quark vs. antiquark (e.g. with charm) would leave very few events. Therefore no simple 
solutions have been found. 

In summary, we have introduced a J2p± measure as a diagnostic of the colour topology 
of hadronic events. The fuzzy nature of hadronic final states somewhat limits the useful- 
ness of the method. In particular, the more drastic effects associated with perturbative 
gluon emission tend to obscure the subtler effects of different colour topologies. Therefore 
the J2p± measure is no panacea, but could still be a useful addition to the (not so large) 
tool box of methods to characterize the nonperturbative stage of hadronic events. Appli- 
cations include three- and four-jet events at Z° energies and, in particular, W pair decay 
to four jets at LEP 2. In the latter process, it could be possible to detect the effects of 
colour reconnection with this approach. Further details on the studies reported here may 
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be found in W£ 
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no reconnection 


with reconnection (model 1) 


algorithm 


1 2 


algorithm 


1 2 


2 


85% 


2 


85% 


3 


68% 65% 


3 


46% 46% 



Table 1: Fraction of agreement in jet pair identification between the three algorithms. 



(a) (b) 

Figure 1: The distribution of J2p.l/N for a: worst— right, and b: maximal— minimal, for 
Lund string (full) and independent (dashed) fragmentation respectively. 
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(a) 



(b) 



(c) (d) 

Figure 2: Distribution of A = (X}p±)worst — (X)p±)best- Dashed lines are reconnected events 
according to models a: 1 (intermediate), b: 2 (bag model), c: 3 (type II superconductor), 
and d: 4 (A-measure), while full lines always are without reconnection. 
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